Gapless points of dimerized quantum spin chains: analytical and 

numerical studies 
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We study the locations of the gapless points which occur for quantum spin chains of finite length 
(with a twisted boundary condition) at particular values of the nearest neighbor dimerization, as 
a function of the spin S and the number of sites. For strong dimerization and large values of 
S, a tunneling calculation reproduces the same results as those obtained from more involved field 
theoretic methods using the non-linear er-model approach. A different analytical calculation of the 
matrix element between the two Neel states gives a set of gapless points; for strong dimerization, 
these differ significantly from the tunneling values. Finally, the exact diagonalization method for a 
finite number of sites yields a set of gapless points which are in good agreement with the Neel state 
calculations for all values of the dimerization, but the agreement with the tunneling values is not 
very good even for large S. This raises questions about possible corrections to the tunneling results. 

PACS numbers: 75.10.Jm, 75.10.Pq, 73.43.Nq 



I. INTRODUCTION 

One-dimensional quantum spin systems have been 
studied extensively for many years, particularly after 
Haldane predicted theoretically that Heisenberg integer 
spin chains should have a gap between the ground state 
and the first excited state [l| , and this was then observed 
experimentally in a spin-1 system Haldane's analy- 
sis used a non-linear er-model (NLSM) which is a field 
theoretic description of the long-distance and low-energy 
modes of the spin system |3|, t& |5l |y, [jj . 

Although the NLSM approach is supposed to be accu- 
rate only for large values of the spin S, it is found to be 
qualitatively correct even for small values of S. For in- 
stance, if there is a dimerization in the nearest neighbor 
Heisenberg couplings, taken to be 1 and k alternately, the 
NLSM predicts that there is a discrete set of values of k 
lying in the range < k < 1 for which the spin chain 
is gapless; these correspond to quantum phase transi- 
tions. Further, the number of gapless points is predicted 
to be the number of integers < S + 1/2; in particular, 
the undimerized chain (with k — 1) is a gapless point 
if S is a half-odd-integer. Numerical analysis shows this 
to be true for values of 5 up to 2 |S IS ES E3; how- 
ever, the numerically obtained values of K at the gapless 
points do not agree well with the NLSM values, ft there- 
fore appears that there must be corrections to the NLSM 
analysis for small values of S. 

The NLSM approach that has been used so far to find 
the gapless points is based on certain properties of a field 
theory in two dimensions (one space and one time); the 
gapless points occur when the coefficient of a topological 
term is given by it modulo 2tt (as will be discussed in 
Sec. II). Although there are arguments justifying this 
criterion E3 j there does not seem to be a simple physical 
picture behind it. One of the aims of our paper will 
be to provide a picture based on tunneling between two 



classical ground states. 

Numerically, there are different ways of finding the 
gapless points for a dimerized spin chain. One of the most 
accurate ways is based on exact diagonalization studies 
of a finite spin chain with a twisted boundar y co ndition, 
to be specified more precisely in Sec. II |lTj,|lj|. In the 
presence of the twist, it is found that the gap between the 
ground state and first excited state vanishes at a value 
of k which is a function of the number of sites 2N (we 
use this notation since the number of sites will always be 
taken to be even). We will be mainly interested in finite 
system sizes in our work; however, it is known from con- 
formal field theory ESESE1ESE1 that the locations 
of gapless points for the infinite system can be found very 
accurately by finding the locations of those points for fi- 
nite systems, and then extrapolating to N — > oo. 

In this paper, we will use three different approaches to 
study finite systems with a twisted boundary condition, 
in order to find the gapless points as a function of S 
and N. The first two approaches will be analytical; they 
will be based on the idea that in the presence of a twist, 
the system has two classical ground states, called Neel 
states, which are degenerate. The degeneracy may be 
broken in quantum mechanics by tunneling. However, if 
the tunneling amplitude is zero, we get a gapless point 
in the sense that the lowest two states have the same 
energy. The third approach will be numerical and will 
be based on exact diagonalization of finite systems using 
various symmetries. 

In Sec. II, we will first define the dimerized quantum 
spin chain and review how it can be described using a 
NLSM field theory. We will then describe the twisted 
boundary condition for a finite system. In Sec. Ill, we 
will describe our first analytical approach. This is based 
on a tunneling calculation for a chain with a finite num- 
ber of sites. For reasons explained below, this method 
is limited to small values of k. We will see that the ex- 
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pressions for the scale of the gap and the locations of the 
gapless points agree with those obtained from the NLSM 
field theory; this is remarkable because our analysis will 
be based only on quantum mechanical tunneling in a fi- 
nite system, while the field theoretic analysis is based on 
a renormalization group equation and the presence of a 
topological term. Our second approach, described in Sec. 
IV, is based on a direct quantum mechanical calculation 
of the matrix element between the two Neel states to low- 
est order in the Hamiltonian; this method works for all 
values of n. We will see that the locations of the gapless 
points obtained by the two methods differ substantially 
for small k. We make some speculations about how the 
tunneling results may be corrected. In Sec. V, we use 
exact diagonalization to find the gaps and the locations 
of the gapless points as functions of K for different val- 
ues of 5 and N. We find that the numerical results for 
the locations of the gapless points agree quite well with 
those found by the second analytical method, and there- 
fore disagree with those found by the tunneling approach. 
In Sec. VI, we will summarize our results. 



II. DIMERIZED QUANTUM SPIN CHAIN 

A. Field theoretic description 

In this subsection, we will briefly review the NLSM 
field theory for a dimcrized spin chain with an infinite 
number of sites. The Hamiltonian is given by 



H 



E 



[5 2 



5 2 



S21 ■ S21+1} 



(1) 



where we have spin 5 at every site, and k > 0. This de- 
scribes an Heisenberg antiferromagnetic spin chain since 
all the couplings are positive. We will set K = 1, so that 
Sf — 5(5 + 1). (We will specify the appropriate bound- 
ary conditions when we discuss finite systems below) . In 
many papers, the nearest neighbor couplings are taken to 
be 1 + 6 and 1 — 5, instead of 1 and n. After a re-scaling 
of the Hamiltonian, we see that the two parameters are 
related as 



1 - 6 
~1+~S 



(2) 



In the classical limit 5 — > 00, the ground state of Eq. 
l(T| is given by a configuration in which all the spins at 
odd sites point in the same direction while all the spins at 
even sites point in the opposite direction. This motivates 
us to define a variable 



S2n—1 — $2n 

25 



(3) 



where x = 2na denotes the spatial coordinate and a is 
the lattice spacing; x becomes a continuous variable in 
the limit a — > 0. In the classical limit, d> becomes a 



unit vector in three dimensions; the model is called the 
NLSM because of this non-linear constraint. One can 
then derive an action in terms of the field 4>(x,i); this 
takes the form 0, |j| 



S 



dtdx 



1 /d0 2 
2cg ( dt ' 



2g [ dx> 



4m 



— / dtdx 6 ■ — x — — 



where 



9 

and 9 



dt 



2aS \[k , 

1 1 + K 

5 



dx 



(4) 



4tt5 
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(5) 



The parameters c, g and 9 denote the spin wave velocity, 
the strength of the interactions between the spin waves 
(even though the first two terms in (Jlj are quadratic, it 
describes an interacting theory because of the non- linear 
constraint on <fi), and the coefficient of a topological term 
respectively. One can show that the term multiplying 9 
in is topological in the sense that its value is always 
an integer. 

It is known that the system governed by Eq. J3J is 
gapless if 9 = it modulo 2tt and g is less than a critical 
value P, 0, ■ This implies that the theory is gapless if 
45k/ (1 + k) = 1, 3, • • • and g is small enough. In partic- 
ular, this means that in the range < k < I, there are 
a discrete set of values of k for which the system is gap- 
less; the number of such values is given by the number 
of integers < 5+1/2. For all values of 9 7^ it modulo 
2ir, the system is gapped. For 9 = modulo 2tt, the 
gap is given by exp(— 2n/g). This follows from the fact 
that the interaction g effectively becomes a function of 
the length scale L and satisfies a renormalization group 
equation of the form dg e ff/d\aL = g2ff/(2Tt). This im- 
plies that g e ff(L) becomes very large at a length scale 
given by Lq ~ a exp(2Tr/g), where g = g e ff(a) is given 
in Eq. This is the correlation length of the system; 
the energy gap is related to the inverse of this length, 
namely, AE ~ cexp(— 2n/g). 



B. Twisted boundary condition 

In this subsection, we will study the same model as in 
Eq. Q but with a finite number of sites going from 1 
to 2N. Although a periodic boundary condition would 
appear to be the simplest, it turns out that a more useful 
boundary condition is one with a twist 

[HIHGl. We 

define the Hamiltonian to be 

N 

H = J]] ['^2n-l^'2n + ^2n-1^2n + ^2n-1^2n] 
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n=l 

+ K [ — ^2N^1 ~ ^2N^1 + S^ffSf] ■ (6) 

Note that the bond going from site 2N to site 1, to be 
called the bond (27V, 1) for short, has a minus sign for 
the xx and yy couplings. We will call this a twisted 
boundary condition; it is equivalent to rotating the x 
and y components of Si by tt about the z axis just for 
that bond. 

An advantage of the twisted boundary condition is that 
classically, the Hamiltonian in © has exactly two ground 
states, namely, (i) S 2 n-i — (0,0,5) and S% n — (0,0,-5) 
for all n, and (ii) S 2n -i = (0, 0, -S) and S 2 „ = (0, 0, S) 
for all n. We will call these two Neel states Ni and 
respectively. (This is in contrast to the case of periodic 
boundary conditions where there is an infinite family of 
classical ground states because S^n-i = — S^n can point 
in any direction). If the classical degeneracy is broken 
quantum mechanically in any way, there will be a gap 
between the lowest two states of the system, while if the 
degeneracy remains unbroken, the system will be gap- 
less. In Sees. Ill and IV, we will describe two ways of 
analytically studying whether the degeneracy is broken. 

Let us now describe the various symmetries of the 
Hamiltonian in Eq. ©. Although the total spin is not a 
good quantum number because of the twist on one bond, 
s tot = J2 n ( S 2n-i + s 2n) is a g ood quantum number. 

Eq. © satisfies the duality property H{n) = kH(1/k), 
where H is related to H by an unitary transformation. 
[The unitary transformation is required because the twist 
only exists at the bond (27V, 1) whose strength is re. By 
applying a rotation Sf — ► —Sf and Sf — > —Sf, we can 
move the twist from the bond (27V, 1) to the bond (1, 2).] 
Since a unitary transformation does not affect the spec- 
trum, we conclude that if there is a gapless point at a 
value re, there must also be a gapless point at 1/re. 

Next, we define the parity transformation P as reflec- 
tion of the system about the bond (2N, 1), namely, 

Si <-> S2N+1-1 ■ (7) 

This is a discrete symmetry of the Hamiltonian H , and all 
eigenstates of H will be eigenstates of P with eigenvalue 
±1. For any value of re, we find that the ground and 
first excited states always have opposite values of the 
parity. We will show below that the relative parity of the 
ground states in the limits k->0 and re — > 00 is given by 
(— 1) 2S . For integer S, this will imply that in the range 
< re < 00, the number of crossings between the ground 
state and the first excited state (and hence the number 
of gapless points) must be even. But for half-odd-integer 
S, the number of such crossings must be odd; combined 
with the duality re — > 1/re, this implies that there must 
be a crossing and therefore a gapless point at re = 1 (this 
is a self-dual point). 



To prove the statement about the relative parities of 
the ground states at re — > and re — > 00 being given by 
(— 1) 2S , we proceed as follows. We first observe that if 
there are only two spins 1 and 2 governed by the un- 
twisted Hamiltonian h u = SfS% + SfS% + Sf S§, then 
under the reflection 1 <-> 2, the ground state has the 
parity (— 1) 2S , while the first excited state has the par- 
ity (— 1) 2S+1 . But for the twisted Hamiltonian fit = 
— Sf 5f — Sf S*2 + £f > tne ground state and first excited 
state have parities equal to 1 and —1 respectively. (This 
can be proved using the Perron-Frobenius theorem for a 
real symmetric matrix). We now consider the entire sys- 
tem with 2Af sites. In the limit re — > 0, the ground state 
of the system is given by a direct product of ground states 
over the bonds (1,2), (3,4), (27V- 1,27V). Since there 
are N dimers, the parity of this state is (— 1) 2SN , while 
the parity of the first excited state is (—1) 2SN+1 . On 
the other hand, in the limit re — * 00, the ground state of 
the system is given by a direct product of ground states 
over the bonds (2,3), (4,5), (2JV,1). Under parity, 
the parity of this state is (_i) 2SAr + 2S ; while the parity 
of the first excited state is (_i) 2 S' Ar + 2S + 1 . A comparison 
between the ground states in the two limits shows that 
they have a relative parity of (— 1) . 

III. TUNNELING APPROACH TO THE FINITE 
SPIN CHAIN 

In this section, we will study the model defined in Eq. 
© using a tunneling approach. We are interested in the 
limit S — > 00 and re — > 0, such that kS is of order 1. We 
will compute the action of the system and use that to 
compute the tunneling amplitude between the two Neel 
states. 

In the limit re — > 0, the system consists of decoupled 
dimers whose energy levels are easy to compute. (We will 
assume in this section that there are at least three dimers, 
i.e, N > 3). For the dimer on the bond (2n — 1, 2n), we 
define the variables 

7* S2n— 1 — S2n 



In — S 2 n-1 + 5*2,! • (8) 

These variables satisfy the relations 

1 P 

& = 1 + - - 

<Pn s 4S 2 > 

4>n ■ L = , 

[c , ft] = * s mn t abc 0« > 

c 

\ja 76] • s abc jc 

c 

and [C , 0»] = * S mn £ e abc l c n . (9) 
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Several simplifications occur in the classical limit S — > 
oo. Firstly, in the Neel state, I is equal to zero while <j> is 
an unit vector; the latter is clear from the first equation 
in J§J) . We will therefore set = 1 exactly. Secondly, we 
will take the commutator [4^, 4> b n ] = due to the fourth 
equation in Finally, given the second and third equa- 
tions in @ , we can obtain the momentum which is canon- 
ically conjugate to 4> n , namely, l n = <f> n x II„, which 
satisfies the commutation relation 



[^m > n n] = i S mn S a b 



(10) 



Using Eq. |JSj, we can write the Hamiltonian in JIjJl in 
terms of <j> and IT, and then obtain the Lagrangian as 



N 



L = £ ^ . 4. - H 



(11) 



We eventually find that 

N 



N-l 



L = y 1 { ^Lf + kS 2 y 

^ 2 K dt ' ^ 



+ KS 2 [ -r N 4>l ~ <t> y N <t> V X + 4>N<I>1 



N-l 



kS \ - , dcj) n 
2 ^ ( dt dt 



n+l 



) • 4>n X <A 



n+l 



kS 



+ (' 



6 w " 



dt dt 

+ + ^)(-«r# + M ] 

+ n 2 S 2 [terms of fourth order in <f> n ] . (12) 

The terms in the third line of this Lagrangian are what 
give rise to the topological term in Eq. Q in the con- 
tinuum limit. In the last line of (|12f> . the terms of fourth 
order in </>„ are chosen in such a way that when we com- 
pute the Hamiltonian from it, it agrees with Eq. JSJ. 
We will now see why these fourth order terms are not 
important. 

In the limit k — > 0, S — > oo and kS of order 1, we 
can scale the time t by a factor of y/S to show that only 
the first two lines of Eq. (|12|) contribute to the Euler- 
Lagrange equations of motion (EOM); this is a major 
simplification. To compute the tunneling amplitude be- 
tween the two Neel states, we will find the solutions of the 
EOM in imaginary time. We then find that the tunnel- 
ing amplitude comes with a phase which arises from the 
third through sixth lines of Eq. I|12|) ; thus these terms are 
important even though they do not directly contribute to 
the EOM. The fourth order terms in the last line of i|12|) 
do not contribute to either the EOM or the phase, and 
we will therefore ignore them henceforth. 



In imaginary time (denoted by the symbol r), the ac- 
tion takes the form 



Si = 



N 



N—l 



^ 2 y dr 

n—l 



<t>n ■ 4>n+l) 



+ kS 2 [ 



<Pn+l 

dr 



) ■ 4>n X <t> n 



+ 1 



. KS d$ n 
n—l 

-i ^ t - + mm 

+ - d -§){-^ - mi) 

+ + ^x-« + 



(13) 



(We have introduced a constant kS 2 N so that the action 
vanishes for each of the two Neel states) . The tunneling 
amplitude will be given by the sum of exp(— Si) along all 
the paths of extremal action which join the Neel states. 
We will now determine these extremal paths. 

Let us use polar angles to write the variables <j> n — 
(sin a n cos @ n , sin a n sin f3 n , cos a n ). The Neel states 1 
and 2 are given by a n — for all n and a n — tt for 
all n respectively. We now solve the EOM following from 
the first two lines of Eq. Ijl3|l i n order to obtain the 
paths going from state 1 to state 2. We will not write 
the EOM explicitly here, but directly present the solu- 
tions. We find that the two paths which have the least 
action are given by 

A: a n (r) = a(r) and f3 n — Po + (nn/N) for all n, and 
B: o„(t) = cc(t) and f3 n = [3q ~ (nir/N) for all n, 
where /3o is an arbitrary angle. In both cases, the function 
a(r) satisfies the EOM 



= nS 2 S in(2a)(l - cos-), (14) 

with the boundary conditions a{— oo) = 0, cc(oo) = tt, 
and da(±oo) / dt = 0. This implies that 



da . . 7r 

— - = 2b V n sma sin — 
dt v 2N 



(15) 



Using this we can evaluate the contribution of the first 
two lines of l|13|) along cither one of the paths joining the 
Neel states. We find that the real part of the action is 
given by 

/OO -| I 

dT [_ (™f + kS 2 sin 2 a (1 - cos ^-)] 
-oo 2 dr N 

= A^SN sin^ . (16) 

We can now evaluate the contribution of the imaginary 
terms in Eq. 1|13JI to the action. We find that for path 
A, they contribute —i2nSNsm(ir/N), while for path B, 
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they contribute i2nSNsm(ir/N). Hence, the total con- 
tribution of the two paths to exp(— Si) is given by 

A - cos(2KSWsm-^) exp(-4\/^SW sin^) , (17) 

up to some pre-factors which are determined by fluctu- 
ations about the classical paths. Since A is the matrix 
element between two classically degenerate states, the 
energy gap between the two states is given by 2|A|. We 
thus see that the gap vanishes if 2kSN s'm(ir/N) is an 
odd multiple of n/2, i.e., if 



4kSN sin — = 7r modulo 2ir 
N 



(18) 



This is the same condition as the one satisfied by the 
parameter 9 in Sec. II A in the limits S, N — ► oo and 
k — ► 0. Further, if AkSN sm{it/N) = modulo 2?r, the 
gap is given by exp(— Re Si) which agrees with the ex- 
pression exp(— 2ir/g) given in Sec. II A for S, N — > oo 
and k — ► 0. Thus, a simple quantum mechanical tunnel- 
ing calculation seems to reproduce the same conditions 
as those obtained earlier by more complex field theoretic 
calculations involving topological terms and a renormal- 
ization group analysis. 

Before ending this section, we should note that there 
are other pairs of paths with extremal action, which have 
the form 

A: a„(r) = a(r) and j3 n = /3q + (pnn/N) for all n, and 
B: a n (r) = a(r) and j3 n = flo — (pnir/N) for all n, 
where p = 3, 5, • • • (going up to the largest odd integer 
< N — 1) labels the different pairs of paths. However, 
the real part of the action of these paths is given by 
4y/nSN sin(p7r/2A), which, for large S, is much larger 
than the expression given in Eq. l|16f) : their contributions 
to the tunneling amplitude are therefore much smaller. 

Finally, we would like to note that it is important 
that the twist in the boundary condition should be by tt, 
and not by any other angle. Even though any non-zero 
twist would lead to two Neel ground states classically, 
the pairs of tunneling paths between those two ground 
states would not have the same real part of the action if 
the twist angle was different from tt. The pairs of paths 
would therefore not cancel each other no matter what 
the imaginary parts of their actions are. (Two complex 
numbers cannot add up to zero, no matter what their 
phases are, if their magnitudes are not equal). 



IV. A SECOND APPROACH TO THE MATRIX 
ELEMENT BETWEEN NEEL STATES 

In the previous section we argued that the gapless 
points of the Hamiltonian in Eq. J§} can be identified 
with the values of k for which the tunneling amplitude 
between the two classical ground states vanish. In this 
section we will calculate this amplitude using an alter- 
nate method. 



The twist on the edge bond (2N, 1) breaks the global 
SU(2) symmetry and thus lifts the continuous degener- 
acy of the classical ground states. With the twist, the 
two degenerate ground states of the Hamiltonian are the 
Neel states which are connected to each other by rotation 
by 7r about the y axis, 

\Nx) = \S,-S,S,---,-S,S,-S) , 
and \N 2 ) = | -S,S,-S,---,S,-S,S) , (19) 

where |{?7ii}) denotes the state with Sf eigenvalue rrij. 
We are interested in the zeroes of the quantity 



T = <JV a |e-^lJVt) 



(20) 



as a function of k. Though the calculation of the above 
matrix element is an interacting many-body problem, 
we can obtain its zeroes 'exactly' in the thermodynamic 
limit. 

First we note that, in the expansion of the exponential, 



-0H 



{-f3H) r 



n=0 



(21) 



the first term which makes a non-zero contribution to T 
has n = 2SN. This is because to take \Ni) to |A 2 ), spins 
belonging to the A-sublattice have to be flipped from \S) 
to | — S), and this requires the action of (S~) for each 
spin. Similarly, the action of (S + ) 2 will take spins in 
the -B-sublattice from 1 — 5) to \S). 

Next, we will calculate the values of k for which 



(N 2 \H 



2SN 



N 1 )=0 



(22) 



Then we will show that as N — > oo, Eq. (1221 implies that 

(A/2|iy 2SJV+fe |iVi) = (23) 

for any finite k. This in turn will imply that T is zero. 

The only term in H 2SN which makes a non-zero con- 
tribution to the left hand side of Eq. 1221) is 



2,5 



We need to count the number of ways in which such a 
term can arise. The contribution comes from terms of 
the following type, 



N 



\2S-m 



(24) 



where < m < 2S. The above term can be obtained in 
(2SN)\ 



(m\) N ((2S - m)\) 



ways and comes with a weight (— l) m K mN . Here we have 
neglected an overall m-independent factor due to the 
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Clebsch- Gordon coefficients arising from the repeated ap- 
plication of S + and S~ operators. The condition in Eq. 
(|2"2"|) then becomes, 



25 N 

]T (-1)™ (mk)) = o 



where a m (n) 



m!(25-m)! 



(25) 



Before proceeding further, we note that the above con- 
dition preserves the duality symmetry under n — > 1/k. 
This is because Eq. (|25|) can be written as, 



„2SN 



(-l) 2S ^(-l) m (a m (l/ K ) 



N 







(26) 



Hence, if n* is a solution, so is Thus we can restrict 

ourselves to the range < k < 1. 

Eq. (|25() determines the roots of a polynomial of order 
2SN , which, in general, cannot be solved analytically. 
But it turns out that we can obtain the roots in the limit 
N — > oo. In this limit, depending on the value of n, one 
particular term in the sum is predominant, and the rest 
of the terms can be neglected compared to this. The 
dominant term is determined by, 



max a r , 



(27) 



As k varies from to 1, to* successively takes values 
m* = 0, 1, • • • , S for integer S and m* = 0, 1, • • • , S + 1/2 
for half-odd-integer S. Noting that neighboring terms in 
the sum have opposite signs, Eq. (|25[> can be satisfied 
only when 



Thus the gapless points are given by 

to + 1 



2S 



(28) 



(29) 



where, m = 0, 1, • • • , S — 1 for integer S and to = 
0, 1, • • • ,S - 1/2 for half-odd-integer S. 

To complete our argument that Eq. (|2~9"|> gives the 
gapless points, we need to show that Eq. I|22|) implies 
Eq. I|23[l . To this end, we first note that the non-zero 
contributions to the matrix element from j{ 2SN + k can be 
obtained from the contributing terms in H 2SN , given in 
Eq. JHJ), by adding terms of the form SfSf +1 , S+S^ +1 
or S^~Sj~ +1 . Now the weight coming from the Clebsch- 
Gordon coefficients depend on the order in which the 
terms appear. But formally one can write that, Eq. Ij23|l 
implies 



5](-l) m 6 m , fe (a m ( K )) Ar = 



(30) 



where 6 m ,fe are finite undetermined constants indepen- 
dent of N. As before, in the large- TV limit, the left hand 



side of Eq. i|30|) can be zero only through the mutual 
cancellation of a pair of neighboring terms, i.e., when 



hn + l 



1/N 



(31) 



As N — > co, this reduces to the condition in Eq. H28J1 . 
In other words, the vanishing of (N2\H 2SN \Ni) is a suffi- 
cient condition for the vanishing of (N 2 \H 2SN+k \Ni) for 
any finite k as N — > oo. 

For half-odd-integer spins, k* = 1 is a solution of Eq. 
(|29|l for to = S — 1/2, but it is not a solution for inte- 
ger spins. This is consistent with Haldane's conjecture 
[1| that the uniform chain is gapless for half-odd-integer 
spins and gapped for integer spins. 

Since the identification of gapless points with the ze- 
roes of the transition amplitudes between the two Neel 
states is essentially a semi-classical approximation, we 
expect the formula given by Eq. I|29|) to get better for 
larger values of S. 

As with the tunneling calculation in Sec III, here also 
one can see that it is crucial to have a twist by tt and not 
any other angle. Let us suppose that the twist angle is 
X- Then the Hamiltonian for the bond between the spins 
at sites 2N and 1 will be 

Then the equivalent of the condition in Eq. (|28|l will read 



e lx a m+ i. 



(32) 



Since a m 's are all real and positive, such a condition can 
be satisfied only when x = 7r, in which case Eq. IMl'l) 
becomes Eq. l|2"5)l. 

Finally, let us compare the gapless values of k given 
in Eq. f2"§|l with those given in Eq. (fT%|l in the limit 
N — > oo, namely, 



1/2 



2S* 



(33) 



where, m = 0, 1, • ■ •. [Since Eq. l|18fl was derived under 
the assumption that S — > oo and kS is of order 1, we 
must restrict to to be much less than S in Eq. I|33|) .] We 
see that for large values of S and m«S, the values of 
n* n in Eqs. (2H1) and J33J) are related by a shift of 1/(45) . 
It would be useful to understand more deeply why this 
is so. Heuristically, this discrepancy can be explained by 
postulating that the phase difference between the actions 
for the two paths discussed in Sec. Ill has an additional 
factor of 7r for some reason. Eq. l|18f) would then change 
to AirnS = modulo 2n for N — > oo; this condition would 
be equivalent to Eq. I|29|l for to <§; S. In a different prob- 
lem (tunneling of a charged particle in two dimensions in 
the presence of a large magnetic field) , it was empirically 
found that an additional factor of ir ap pears due to fluc- 
tuations about the tunneling paths [18j . It may be worth 
studying if something similar happens in our problem. 
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V. NUMERICAL RESULTS FOR FINITE 
SYSTEMS 



In this section, we numerically determine the values 
of k for which the Hamiltonian becomes gapless using 
exact diagonalization of finite systems. These results, 
being a direct calculation of the gapless points, will give 
us information about the physical regimes in which the 
analytical methods outlined in Sees. Ill and IV are valid. 

We begin with a brief outline of the method used to 
find the gapless points. As we vary k in Eq. {T|l. it is 
known that the gapless points separate various phases 
whose ground states are represented approximately by 
different valence bond states (see Refs. [n], US an d ref- 
erences therein) . Specifically for a given spin 5, there are 
25 + 1 different phases separated by the 25 gapless points 
for k between and oo. The lowest two eigenstates of an 
untwisted Hamiltonian never cross each other in energy 
even at the transition from one phase to the other; also, 
the ground state always has the same eigenvalue (—1) 2SN 
for the parity P. It is here that we make use of the twisted 
boundary condition in Eq. ©. It has been shown [lfl ] 
that the lowest two eigenstates of this Hamiltonian (both 
of which lie in the 5 t z ot = sector) have different parity 
eigenvalues, and they cross at certain points which, in the 
limit ./V — > oo, become the gapless points of the Hamil- 
tonian in Eq. This means that one can locate the 
gapless points of the Hamiltonian in Eq. Q by studying 
the crossing of the two lowest eigenvalues of the Hamil- 
tonian in Eq. © in different parity sectors. This enables 
us to find the gapless points without having to consider 
two very closely spaced eigenvalues lying within a single 
symmetry sector. 

We study the Hamiltonian in Eq. © numerically. Be- 
cause of the duality between k and 1/k we restrict our 
studies to < k < 1. We use the Lanczos algorithm to 
diagonalize finite systems from N = 3 to N — 5. Spins 
from 5 — 2 to 3 are studied for TV = 3 to N = 5, but for 
5 = 3.5 to 7 we restrict ourselves to N = 3 and N = 4 
due to computational limitations. 

In Fig. 1, we present two representative data plots, 
for N — 3. The upper plot is for 5 = 3.5 and the lower 
one is for 5 = 4. Plotted on the y axis in both figures 
is E(P = —1) — E(P = +1), the energy difference be- 
tween the two lowest energy eigenstates in the two parity 
sectors. Wherever the plot crosses the x axis we have a 
gapless point of the Hamiltonian in Eq. QJ. The values 
of the crossing points are indicated in the plots. From the 
discussion in Sec. II, we know that the ground state for 
N = 3, 5 = 3.5 should have a parity —1, and for N = 3, 
5 = 4 should have a parity +1. This is indeed borne out 
by the plots in the figure. Moreover, the gapless point 
at k = 1 is also present for 5 = 3.5 as expected. The 
nature of these plots for other spins and different lattice 
sizes is similar. Since our emphasis in this work is on the 
locations of the gapless points, we now turn to analyzing 
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FIG. 1: Energy difference between the two lowest energies in 
the P — — 1 and P = +1 sectors as a function of k, for N = 3. 
The upper figure is for S = 3.5, and the lower figure is for 
5 = 4. The locations of the gapless points are shown. 



those points more closely. [Incidentally, we observe in 
Fig. 1 that the envelope of the magnitude of the gap is 
rapidly decreasing with increasing k; this is in accordance 
with the exponential factor for the tunneling amplitude 
inEq. JTZJ.] 

Figure 2 shows a comparison of the different methods 
used to calculate the gapless points. Plotted on the y axis 
is K , the gapless point closest to the origin for various 
values of the spin, for N = 3. The topmost plot (marked 
by dots) is for values obtained from the numerical calcu- 
lations, called Kn Um - The next plot (crosses) is for values 
obtained from the Neel state calculation in Eq. 125fl . 
K Neei ■ The plot at the bottom (squares) is for values ob- 
tained from the tunneling expression in Eq. (JTSJ, k\ u11 . 
Clearly, the values of k 1 obtained using the Neel state 
calculation are in closer agreement with the numerical 
values than the tunneling values. The tunneling values 
do not converge with increasing 5, even though one is 
looking at data for spin values as large as 5 = 7. On the 
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FIG. 2: The location re 1 of the gapless point closest to zero, 
as a function of the spin, for TV = 3. The results from nu- 
merical calculations (dots), Neel state calculations (crosses) 
and tunneling calculations (squares) are shown. The joining 
lines are guides for the eye. The inset shows the percentage 
variation of the Neel state calculations when compared with 
the numerical results. 
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FIG. 3: The location re* of the gapless point closest to (but 
less than) 1, as a function of the spin, for TV = 3. The results 
from numerical calculations (dots) and Neel state calculations 
(crosses) are shown. The top and bottom parts of the figure 
are for integer and half-odd-integer values of the spin respec- 
tivly. The joining lines are guides for the eye. The insets 
shows the percentage variation of the Neel state calculations 
when compared with the numerical results. 

other hand, the plot of values using the Neel state calcu- 
lations converges much faster. This is made clearer in the 
plot by the inset where on the y axis we have plotted the 
percentage deviation, ((«;* um - K N eel)/ K ium) X 100, of 
the values of K} Neel from the numerically obtained values. 
Figure 3 shows a comparison between the numerical re- 



sults (dots) and those obtained from the Neel state calcu- 
lation (crosses) for re*, the gapless point closest to k = 1. 
[Unlike Fig. 2, we have not shown the tunneling results 
based on Eq. (|18fl ) because that formula for the gapless 
points is not valid when re is close to 1.] The data sets for 
integer and half-odd-integer spins have been plotted sep- 
arately. The plot and the inset at the top are for integer 
spins, and the ones at the bottom are for half-odd-integer 
spins. As before, we see that the agreement with the nu- 
merical results improves as we go to larger spins. We also 
see from the insets that for a given spin, the agreement 
is much better near re = 1 than it was near re = 0. This 
suggests that the Neel state calculation gets better as we 
increase re from to 1. This is something which is seen 
very clearly in the following tables. 





KNeel 


% deviation 


r^tun 


% deviation 


0.083 


0.080 


3.6% 


0.047 


45.4% 


0.190 


0.180 


5.2% 


0.140 


26.3% 


0.306 


0.293 


4.2% 


0.232 


24.2% 


0.438 


0.424 


3.19% 


0.326 


25.6% 


0.773 


0.765 


1.03% 


0.512 


33.8% 


1 


1 


0% 







^nura 


KNeel 


% deviation 




% deviation 


0.077 


0.074 


3.89% 


0.043 


44.2% 


0.174 


0.166 


4.02% 


0.130 


25.3% 


0.281 


0.269 


4.98% 


0.252 


23.1% 


0.400 


0.387 


3.25% 


0.302 


24.5% 


0.536 


0.525 


2.05% 


0.389 


27.4% 


0.695 


0.685 


1.29% 


0.475 


31.7% 


0.887 


0.884 


0.34% 


0.561 


36.8% 



Table 1. Comparison of the values of all the gapless 
points obtained using the three methods of calculating 
them, for S — 6.5 (top) and 7 (bottom). The data pre- 
sented is for N = 3. 

Table 1 shows how the numerical results, the Neel state 
calculations and the tunneling results compare for all the 
gapless points. The table shows the values of re at which 
the Hamiltonian in Eq. Q is gapless for S — 6.5 (top) 
and S = 7 (bottom), for TV = 3. The percentage devia- 
tions as defined earlier are also shown for the Neel state 
and tunneling calculations. As conjectured after looking 
at Fig. 3, we see that the Neel state calculation gives 
successively better approximations to the actual gapless 
points as we go further from the origin re = 0. We have 
shown the tunneling values for all the gapless points only 
for completion; the formula given by Eq. (|18fl is valid 
only for nS of order 1. We again see that these values 
have much larger percentage deviations from the values 
obtained from the other two methods, even though the 
values of the spins considered are quite large. 

We now look at how the values of re at the gapless 
points change with N and see how the N — > oo values 
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1/(2N) Z 

FIG. 4: Variation of re 1 with N for S = 2.5, 3, 3.5, and 
4. The numbers at the left of each graph are the N — > oo 
extrapolated values. 

compare with those given by Eq. H29[). We take re 1 as an 
example. Figure 3 shows the behavior for S — 2.5, 3, 3.5, 
and 4. We find the N — > oo values by extrapolating the 
best fits obtained by fittin g th e data to even polynomials 
in l/(2iV) 2 following Ref. jlj. All the data are for N = 
2,3, 4 and 5 (for 5 = 2.5 and 3). 

The extrapolated values of re 1 in the N — > oo limit 
and the corresponding values obtained from Eq. i|29|) 
(in parentheses) for S — 2.5, 3.0, 3.5 and 4.0 are 
given by 0.223 (0.200), 0.179 (0.167), 0.146 (0.143) and 
0.124 (0.125) respectively. Clearly, the agreement be- 
tween the two sets of values gets better for larger values 
of the spin. 

VI. CONCLUSIONS 

We have used three different techniques to find the 
gapless points of a dimerized spin-.!? chain with a finite 
number of sites and with a twisted boundary condition. 
The first technique uses a tunneling approach which is 
expected to be valid in the limit S — > oo, re — > and kS of 
order 1. Remarkably, we find that a quantum mechanical 
tunneling calculation reproduces the same expressions for 
the locations of the gapless points and the gap as those 
obtained by more involved field theoretic techniques. 

However, a direct numerical study of the gapless points 
shows a systematic deviation from the tunneling results 
in the limit discussed above. It would be useful to know 
why the tunneling results differ systematically from the 
numerical results in this limit. One possible idea is to 
examine if an additional factor of n appears in the fluctu- 
ation pre-factor of the tunneling amplitude as mentioned 
at the end of Sec. IV. 

In view of the discrepancy between the tunneling and 



numerical results, we have presented a second analyti- 
cal derivation of the gapless points which is based on a 
calculation of the matrix element between the two Neel 
states to lowest order in powers of the Hamiltonian; this 
derivation is expected to become more accurate as the 
number of sites becomes large. We find that the results 
obtained by this approach agree much better with the 
numerical results than the tunneling results, even in the 
limit re — ► 0. It may be instructive to understand in more 
detail why there is such a good agreement between this 
relatively simple analytical calculation and the numerical 
results. 

One of the features of the numerical results shown 
above is that the Neel state calculation always underesti- 
mates the values of re which correspond to gapless points, 
while all the time getting closer to the actual values with 
increasing S (given N) or increasing N (given S). This 
may indicate positive corrections of order 1/N and 1/5 
to the formula obtained from the Neel state calculation. 
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